No-go theorem in many body dissipative particle dynamics 
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Many body dissipative particle dynamics (MDPD) is a particle-based simulation method in which 
the interaction potential is a sum of self energies depending on locally-sampled density variables. 
This functional form gives rise to density-dependent pairwise forces, however not all such force laws 
are derivable from a potential and the integrability condition for this to be the case provides a 
strong constraint. A strategy to assess the implications of this constraint is illustrated here by the 
derivation of a useful no-go theorem for multicomponent MDPD. 



PACS numbers: 61.20.Ja, 05.20.Jj 



Dissipative particle dynamics (DPD) has attracted 
a lot of interest in its possibilities for modelling soft 
condensed matter Q, Q. It is characterised by pair- 
wise soft repulsive forces with a pairwise momentum- 
conserving thermostat Q. Conceived somewhat later, 
many-body dissipative particle dynamics (MDPD) holds 
much promise as a second generation method [4l4ll|. al- 
though it has mainly been applied to vapour-liquid coex- 
istence and free surface simulations [7H11J] . In the present 
note I explore the consequences for MDPD of allowing 
an arbitrary local density dependence into the pairwise 
repulsive forces. In particular if one requires that a po- 
tential exists, so that the forces are conservative, the al- 
lowed functional form of the pairwise forces is severely 
constrained. I will here outline a strategy for assessing 
the implications of this constraint, illustrated by the gen- 
eration of a no-go theorem for a multicomponent MDPD 
force law that has been suggested in the literature 0, [l(| ■ 
The present result also demonstrates that it is by far bet- 
ter to proceed from the potential to the forces, rather 
than the other way around. 

Leaving aside the DPD thermostat, which has been 
well described elsewhere [L2} . standard DPD is charac- 
terised by pairwise repulsive forces of the form 



Fij = Awc(nj)r 



n ) * n 



(1) 



where Fy is the force acting between the zth and jth par- 
ticles (at positions r^ and r^), A is a repulsion amplitude, 
wc( r ) is a weight function, Tij — \vj — r;| is the spatial 
separation, and f^ = (rj — Vi)jrij is a unit vector along 
the line of centres. I shall assume the weight function 
has compact support (i. e. wc — for r > r c where r c is 
a cut-off distance) but for the present purposes it is not 
necessary to specify the exact functional form. 

In contrast, MDPD starts from a rather different view- 
point. The potential energy in^DPD is a sum of density 
dependent one-body terms 



tf({*<}) = Et«(ft), 



(2) 



where the one-body terms depend on local densities, 

Pi = T,ijtj^ P (nj)- (3) 

From the potential one can derive the force law, 
dU 



Or, 



= E 



ijtj F V > 



where 



-[u'ip^ + u'ip^w'pin^tij 



(4) 



(5) 



Here w p (r) > is another weight function, also with 
compact support. Since the potential in Eq. ([2]) is a reg- 
ular function of the particle positions, MDPD avoids is- 
sues that otherwise plague density-dependent pair inter- 
actions [13| , although if the forces are not purely repulsive 
one should take care to ensure thermodynamic stability 
according to the criteria devised by Ruelle fl4j |. 

It is clear that the choice u(p) = A~p/2 and w' (r) = 
—wc(r) brings Eq. ([5]) into agreement with Eq. ([1}. 
Hence standard DPD is just a special case of MDPD. 
Note that this may imply Jd 3 rw p (r) ^ 1 but abandon- 
ing this normalisation requirement leads to a considerable 
notational simplification by eliminating unnecessary pref- 
actors. Another example is u(p) = B~p 2 /2. Again setting 
wc = — w'ni this generates the force law 



B{Pi+ p 1 )w c {r lJ )r l . 



(6) 



This force law (with B > 0) in combination with the 
standard DPD force law of Eq. ((T|) (with A < and a 
larger cut-off) has been extensively used for free surface 
simulations. For a recent review see Ghoufi et al. [ill ] 

Frequently DPD is applied to multicomponent systems 
and Eq. ^ is generalised to 



F^ = AijWc(nj)r. 



ij j l ij 



(7) 



where A.y is a matrix of repulsion amplitudes. It is nat- 
ural to consider a similar generalisation of Eq. (|6]), 



B ij(Pi + Pj)wc(rij) 



y 1 x *j • 
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(8) 

This has been proposed in the published literature |9l.[l0|. 
but my claim is that such a force law is not conservative 



unless Bij is a constant matrix. This is the no-go theorem Thus we see the Dij vanish if and only if Bn = B\ 3 
of the title. (In fairness to the authors of Refs. |9|, Il0| . 
they actually only use the By = B case.) 

How can the no-go theorem be proved? Hopefully it is 
obvious that a sufficient condition for a many body force 
law to be conservative is to display an explicit potential. 
For the present problem, if By = B, such a potential 
is provided by Eqs. © and (O with u{p) = B~p 2 /2 as 
already stated. A necessary condition for a force law to 
be conservative is that the 'Maxwell relation' 



B23, since by choice all the w' (xy) factors are strictly 
negative. Moreover we can pick any three particles for 
this argument. Hence the Maxwell relation fails in at 
least in a subset of configurations, unless By is a constant 
matrix. This completes the desired proof of necessity. 



<9r, 






(9) 



is satisfied. This follows from the first part of Eq. (QJ. 
The application of this to a general configuration of N 
particles is not straightforward since both sides of Eq. (|9]) 
contain multiple sums. However Eq. Q should hold for 
any configuration of particles, so we can choose a con- 
figuration at our convenience. For the present problem 
therefore, let us select one which contains an isolated 
collinear triplet of particles. Without loss of general- 
ity we can place the particles on the x-axis at positions 



X\ < X2 < x 3 . I shall define Xy = Xj 



We can fur- 



ther suppose X13 < r c so that all particles interact. The 
pairwise forces are (setting wc — -iul) 



(10) 



The local densities are ~p 1 = w p (xi2) + w p (xi3), p 2 — 
Wp{xi2) + Wp{x2s,) 7 p 3 = w p (xi 3 )+w p (x23); and the total 
forces are F x = F 12 +F 13 , F 2 = F 2 i+F 23 , F 3 = F 31 +F 32 . 
Let us define Dy = dFi/dxj — dFj/dxi. By explicit 
calculation I find D12 = D 23 = D31 where 



F12 = 


—F21 — 


-£12(^+^2)^(2:12) 


Fi 3 = 


—F 3 \ = 


-Bi 3 {-p 1 + p 3 )w' p (x 13 ) 


F23, — 


— F 3 2 = 


-B 2 3{P2+P3)'W' P {X2 3 ) 



F>i2 = (Bi 3 - B 23 ) w' p {x 13 ) w' p (x 23 ) 

+ (B12 - B 13 ) w' p (x 12 ) w' p (x 13 ) 

+ (B12 - B 23 ) w' p (x 12 ) w' p {x 23 ) 



(11) 



An analogous argument does not go through for Eq. ([7]) 
since in that case we can exhibit an actual potential, 
namely U = J2i>j Ayw p (ry), although this does not re- 



duce to a sum of self energies unless Ai 



Aj. 



Of course one can additively combine force laws, so it 
is possible to have conservative multicomponent MDPD 
based on Eqs. (|6]) and ([7]). Other approaches have been 
described by Trofimov et al. [5J and Merabia et al. [8|. 



Why is it so important that the force law be conserva- 
tive? The answer is that for many applications one wishes 
to use the machinery of equilibrium statistical mechanics 
and thermodynamics |15| . which requires the existence 
of a potential so that the steady state is given by Gibbs- 
Boltzmann. Of course the absence of a potential does not 
preclude the existence of a non-equilibrium steady state. 
Hence the effects of a non-conservative force law could be 
quite subtle, a bit like a failure to satisfy detailed balance 
in a Monte-Carlo simulation [l|. 



Here the general strategy to prove the no-go theorem 
has been to identify a convenient particle configuration 
for which the sums implicit in Eq. (J9j> become manage- 
able. Presumably this could be applied in other cases 
too, but it would seem that at least N > 3 particles are 
required since a central force between a pair of particles 
can always be integrated to a pair potential. The con- 
verse implication is that a similarly defined potential for 
assembling N > 3 particles would depend on the assem- 
bly path, unless Eq. © holds. 
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